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Abstract 



-The molecular ion HfF + is the chosen species for a JILA experiment to measure the electron electric dipole moment 
[(eEDM). Detailed knowledge of the spectrum of HfF is crucial to prepare HfF + in a state suitable for performing an 
-eEDM measurement [1]. We investigated the near- infrared electronic spectrum of HfF using laser- induced fluorescence 
1> (LIF) of a supersonic molecular beam. We discovered eight unreported bands, and assign each of them unambiguously, 
-four to vibrational bands belonging to the transition [13.8]0.5 A1.5, and four to vibrational bands belonging to the 
[transition [14.2]1.5 <— A1.5. Additionally, we report an improved measurement of vibrational spacing of the ground 
state, as well as anharmonicity u e x e . 
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q 1. Introduction 
(— > . 

03 HfF+ is one of the candidate species for experiments to 
& measure the electron electric dipole moment (eEDM). In 



particular, the long-lived mctastable Ai state of the ion 



^[possesses both a long coherence time [2, 3] and a strong 
pH^effective electric field [2, 4], offering enhanced sensitivity 
O |to an eEDM signal of perhaps as many as two orders of 
'— '-magnitude [1] over the current experimental limit [5]. The 
[degree to which an electron possesses an electric dipole 
-moment is deeply connected to the phenomenon of charge- 
parity violation in particle physics. A more sensitive mea- 
-surement of the eEDM will stand to constrain the param- 
etcrs of various so-called Beyond Standard Model theories 
-of high energy physics. 

A crucial aspect of the eEDM measurement is prepar- 
ing a sample of HfF + in the state with the highest eEDM 
^ sensitivity, in a given magnetic sublevel of the J = 1 man- 
• • -ifold of 3 Ai v = 0. The proposed mechanism for this is 
I two-color autoionization directly into the science state [6]. 
^ -To this end, we require detailed knowledge of the spectrum 
S-h [of HfF to identify intermediate states for the autoioniza- 



o 
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tion process. 

Spectroscopy of HfF has previously been done by Adam 
et al.[7]. In addition to the present work, we have done 
spectroscopy on HfF using resonantly-enhanced multi-photon 
ionization [8]. Also of relevance is the work done by Barker 
et al.[3], which is the first spectroscopy of the low energy 
states of HfF + , as well as other work by our group inves- 
tigating the higher energy states of HfF + [9, 10]. 
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2. Experiment 

We use a supersonic beam apparatus similar to that 
which is described in Ref. 6. We hold a mixture of 99%Ar 
+ 1%SF 6 at 100 psi backing pressure. A pulsed valve opens 
for 140 /is and the gas undergoes expansion, achieving su- 
personic velocities after entering the ablation chamber. A 
30 mJ pulse of focused light at 1064 nm from an Nd:YAG 
laser ablates material from a Hf target rod 0.14 in. in 
diameter. The ablation plume is entrained in the super- 
sonic beam and chemically reacts to produce HfF, HfF + , 
and other products. The molecular beam passes through a 
beam skimmer 11 cm from the point of ablation, and then 
passes through a second skimmer 9.6 cm further down- 
stream. The diameters of the first and second skimmers 
are 3 mm and 2 mm respectively, and the region between 
the skimmers is differentially pumped. At a distance of 8.5 
cm downstream from the second skimmer the molecules 
are interrogated using laser- induced fluorescence (LIF) . A 
130 mJ pulse of light, with a linewidth of 0.04 cm -1 , from 
a tunable dye laser excites the molecules, which then emit 
fluorescence which is collected with a gold-coated spherical 
mirror and focused with an off-axis parabolic mirror onto 
an R3896 Hamamatsu photomultiplicr tube (PMT). The 
fluorescence lifetime is long compared to the transit time 
of 15 its. In order to temporally gate the prompt scattered 
light of the fluorescence laser, there is a 6 /lis delay between 
the laser firing and the beginning of the 9 its window where 
the PMT is active. The wavelength range of our scans was 
covered by a single dye, LDS722. 

For the survey spectroscopy in Fig. 1, the dye laser 
is calibrated to Argon lines in an optogalvanic cell, and 
the laser frequency is determined by interpolating the laser 
grating position between measurements made with a waveme- 
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Figure 1: Low resolution survey spectroscopy, with the positions of the eight vibrational bands reported in Table 1 and the Hf 3 F2 — > 3 Di 
transition labeled. 



ter at the beginning and end of the scan. For the high- 
resolution spectroscopy of the individual bands, for ex- 
ample see Fig. 2, the laser frequency is measured at each 
point in the scan with two separate wavemeters, one with a 
high absolute accuracy and the other with high resolution, 
specified to be 10 MHz. The high-resolution wavemeter is 
calibrated daily to an external-cavity diode laser locked to 
a 87 Rb transition at 384.227982 THz. 



3. Results and discussion 
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Figure 2: The [14.2]1.5 - XI. 5(0, 0) band of HfF with ground and 
excited state fit to the Hamiltonian in Eq. 1. This demonstrates a 
Av = transition with an Q' = 3/2 upper state. The fit is inverted 
and displayed under the data. 



analysis we looked for and subsequently found an eighth 
band outside of the original survey range. 

In the eight high-resolution scans we see four bands 
with unresolved isotope splitting, consistent with Av = 0, 
and visibly distinct P, Q, and R branches. Of these bands, 
on the basis of the presence or absence of certain low J 
lines, two are clearly (ft' = 3/2) <— (il" = 3/2), and two 
are clearly (J7' = 1/2) <— (SI" = 3/2) with considerable 
Si-doubling. Based on values of the rotational constants, 
band origin v, and isotope shifts, we are able to assign 
these four bands and the remaining four less well-resolved 
bands to a common lower electronic state, with rotational 
constants in good agreement with the XI. 5 state identified 
in Ref. 7. The upper state of these bands we assign to 
one of two newly identified electronic states: a state with 
Q' = 3/2 at 14223 cm" 1 , which we label [14.2]1.5, and 
a state with fi' = 1/2 at 13833 cm -1 , which we label 
[13.8]0.5. 

On each of the bands acquired at high resolution we 
perform a non-linear least squares fit to a vibrational band 
contour. The doubling for all O = 3/2 levels is unresolved. 
For the lower states of all eight bands and for the upper 
states of the bands at 13496, 13554, 14159, and 14223 
cm -1 we fit using the following Hamiltonian: 



F(J) = T + BJ(J +1). 



(1) 



The bands at 13760, 13833, 14352, and 14429 cm^ 1 we fit 
to an upper state using the Hamiltonian[ll], 



F a (N) =T + BN(N + 1)-±(N+ 1), 
F b (N) =T + BN(N + 1) + 1 N, 



(2) 



From the survey scan we chose seven vibrational bands 
to acquire at higher resolution, and following preliminary 
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Figure 3: The [13.8]0.5 - XI. 5(0, 0) band of HfF with excited state 
fit to the Hamiltonian in Eq. 2. This spectrum is an example of a 
Av = transition with an Q' = 1/2 upper state with f!-doubling. 
The fit is inverted and displayed under the data. 



with selection rules for the P, Q, and R branches: 



Pa : 


(n' 


4) 


-J" 


= -1 


Pb : 


(n' 




-J" 


= -1 


Qa : 


(n 1 




-J" 


= 0, 


Qb : 


(n' 




-J" 


= 0, 


Ra '■ 


(n' 


-I) 


-J" 


= +1 


Rb : 


(n' 




-J" 


= +1 



(3) 



The absolute parity of the excited state could not be deter- 
mined, so following the example in Ref. 7, we label related 
pairs as a/6 instead of e/ f as determined by parity. We do 
not observe a sufficient number of rotational lines within 
each band to be able to fit a non-zero value to the higher 
order rotational constant D, as such it is not included in 
our fitting models. 

There are five isotopologs of HfF, labeled ™Hf 19 F, (with 
n = 180, 179, 178, 177, and 176) each with corresponding 
reduced mass and respective relative abundances 35.1%, 
13.6%, 27.3%, 18.6%, and 5.26%[12]. For Av = transi- 
tions (for instance, as in Figs. 2 and 3) the isotope struc- 
ture is not resolved in our spectra as it is less than the 
0.04 cm -1 linewidth of the LIF laser. For Aw = +1 and 
Av = — 1 transitions the isotope spacing is in principle re- 
solvable, but with the exception of notably sparse regions 
of the spectrum (Fig. 4) we are unable to assign indi- 
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Figure 4: A relatively uncluttered portion of the [14.2] 1.5— XI. 5(1, 2) 
band highlighting the R(1.5) line, where the 3 most abundant iso- 
topologs of HfF are distinctly visible. Plotted for comparison are 
the predicted positions of these lines calculated from the reduced 
mass of each species and the relative intensities expected from nat- 
ural abundance of each Hf isotope. For low J lines, the intensity of 
fluorescence peaks from the 177 Hf 19 F and 179 Hf 19 F isotopologs may 
be diluted due to hafnium hyperfine structure. 
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Figure 5: The [13.8]0.5 - XI. 5(2, 1) band of HfF with excited state 
fit to the Hamiltonian in Eq. 2. The fit is inverted and displayed 
under the data. For bands of lower quality such as this, our fit 
only includes *V, r;, overall intensity and rotational temperature as 
adjustable parameters. Rotational constants are determined from 
fits to other bands. 



vidual rotational lines due to the overlapping spectrum of 
each isotopolog (as in Fig. 5). 

In light of this difficulty we adopt the following fitting 
strategy. We assume that the rotational constants of each 
isotopolog scale reliably with reduced mass n fi: 



(4) 



where n = 178.494 labels the abundance-weighted mass 
average of Hf, and the quantities n \i = 17.171 amu and 
n B are the reduced mass and rotational constant of this 
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Table 1: Molecular constants of vibrational bands of HfF. derived from global non-linear least squares fits. The errors determined by 
bootstrapping the residuals of the non-linear fits are reported in parentheses as one standard deviation of the last digit. Values reported 
without uncertainties are held constant in their respective fitting routine. As described in the text, r\ is the per-amu isotope splitting parameter, 
and n v is the band origin of fictitious isotopolog n Hf 19 F. The fit value for n B is n B( n [i/ n fj,), and the fit value for n u is *V + r](n — n) 



fictitious mass-avcraged isotopolog. Similarly, we take the 
band origins of each isotopolog to be distributed about 
their averaged value as follows: 



= n v + 77(71 - n), 



(5) 



where 77 is the per-amu isotope splitting parameter. The 
reduced mass of any isotopolog differs from the average 
mass by a fractional amount of no more than 0.2%, jus- 
tifying this linear expansion at the level of accuracy we 
aspire to. To the extent the isotope shift is dominated by 
the vibrational contribution, we expect 77 to be 



'/' 



AV, 



(0) 



e=0 



where AV is the difference in vibrational energy between 
the upper and lower states. It is assumed that at our level 
of accuracy we are not able to observe the isotope shift in 
the spin-doubling parameter 7. 

We use different fitting procedures for Av = bands 
compared to Av = ±1 bands. For Av = bands we 
assume 77 = and fit values of n v, n B' , n B" , tempera- 
ture, instrumental resolution, and for the f2' = 1/2 bands 
also 7. When the same vibrational levels are involved in 
bands with Av ^ 0, we use previously determined values 
of n B', n B" , and 7, and fit a contour in order to extract 
n v and 77 only. The exception to this procedure was for 
the v = 2 levels of the XI. 5 and [13.8]0.5 states. For 
these states we had no alternative but to extract the cor- 
responding values of B from the less well resolved bands. 
We calculate tentative, purely statistical, one-sigma er- 
ror estimates using a bootstrap method on the resampled 
residuals[13]. This method assumes the fitting model is 
perfect, and ignores the possibility of, e.g., nonthermal ro- 
tational populations or other deviations from Honl-London 
intensity predictions. We treat these estimates as indica- 
tive of the relative sizes of the various error bars. From in- 
consistencies observed in values of parameters determined 
redundantly by fits to different bands, we find that a global 
doubling of each error bar is appropriate. A summary of 
the combined results and rescaled one-sigma error bars is 
presented in Table 1. 



Even though the quality of our spectra allows for lit- 
tle redundancy in the determination of the rotational con- 
stants £?, our assignment of vibrational levels is quite solid. 
Observed isotope shifts are consistent with the assigned 
values of Av, and the redundant pair of measurements 
for the B" for XI. 5 v = is consistent to a difference 
of < 0.001 cm" 1 . The two bands with O' = 3/2 and 
v = permit a precise determination of AG"^ (668.641 (7) 
cm -1 ) for the XI. 5 state, and the two bands with O' = 1/2 
and v = 1 give a redundant, if less precise determination 
of the same value. The combined result is in good agree- 
ment with, but much more precise than, the previous result 
670(13) cm -1 reported by Adam et al[7]. 
We also observe one transition from the v = 2 level of 
the XI. 5 state. From this we are able to report a value 
AG3/2 = 663.082(9) cm^ 1 . This allows us to extract the 
anharmonicities Lo e x e of the XI. 5 and the [13.8]0.5 states, 
which are 2.78(1) cm -1 and 2.6(1) cm -1 , respectively. Our 
measured value of the ground state ui e x e is consistent with 
the calculation based on the functional form of a Morse 
oscillator [14], uj1/AD e ~ 2.0 cm -1 , calculated from the 
measured bond-dissociation energy D e = 54000 cm _1 [15]. 
A summary of the molecular constants for different vibra- 
tional levels of the states XI. 5, [13.8]0.5, and [14.2]1.5 is 
in Table 2. 

We are not aware of published theoretical predictions 
for HfF, but the authors of Rcf. 16 provide predictions 
for the isoelectronic species HfCl. They predict three clus- 
ters of states: (i) three low-lying doublet states including 
a ground state 2 A 3 / 2 , (ii) a collection of quartet states 
around 9000 cm" 1 , and (iii) a collection of doublet states 
around 18000 cm -1 . In the experimental work on HfF, 
the authors of Ref. 7, having worked from the assumption 
that HfF structure will be similar to HfCl, speculated that 
the seven doublet levels they observed from 17000 cm" 1 
to 24000 cm" 1 arc part of cluster (iii). If this assumption 
is true, it is plausible that the [14.2] 1.5 and [13.8]0.5 levels 
we have discussed are derived from the nominally quartet 
levels of group (ii) . 
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0.2690(20) 
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14222.928(6) 
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Table 2: Molecular constants for the observed vibrational states of HfF. Values indicated by a dash are where no 7 parameter is included 
in the Hamiltonian describing that state. Values noted by an asterisk are for states where the data is not well enough resolved to obtain a 
precise value of 7 from a fit. At this level of accuracy, 7 is not likely to depend on the vibrational level, and thus we present the measured 
value of 7 for the [13.8]0.5 V = state as a suggested value for the v = 1 and v = 2 levels of the same state. 



4. Conclusion 

We extend observations of the spectrum of HfF using 
LIF on a supersonic beam. After recording a survey spec- 
trum from 13500 to 14500 cm -1 we acquired eight vibra- 
tional bands at higher resolution. We perform rotational 
analysis on these bands, assigning them as vibrational sub- 
bands of two new electronic states of HfF, [13. 8] 0.5 and 
[14. 2] 1.5. From the analysis we are also able to extract 
rotational constants for the electronic states, as well as 
AG" ^ 2 an d w e % e for the ground state. 
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